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function  of  Sn/n^  (fractional  density  depletion),  depending  on  the  shape  of  the  bubble.  Con- 
sequently, high  altitudes  and/or  large  percentage  plasma  depletions  yield  large  vertical  bubble  rise 
velocities.  The  models  are  extended  to  include  arbitrary  electric  fields  in  addition  to  the  gravi- 
tational force. 
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MORPHOLOGICAL  STUDIES  OF  RISING  EQUATORIAL  SPREAD  F BUBBLES 


■ 


I.  INTRODUCTION 

Recently,  experimental  (Kel ley  et  al . , 1976;  Woodman  and  La  Hoz, 
1976;  and  McCl ure  et  al . , 1977)  and  numerical  simulation  studies 
(Scannapieco  and  Ossakow,  1976)  of  the  nighttime  equatorial  F region  have 
given  evidence  for  rising  bubbles  (plasma  density  depletions).  The  ver- 
tical speed  of  these  bubbles  has  been  observed  (McClure  et  al . , 1977)  to 
range  from  very  low  (several  meters/sec  and  even  zero)  to  very  high 
('vlOO  m/sec  or  more)  velocities.  Hudson  [1977]  has  attempted  to 
describe  these  rising  bubbles  within  the  context  of  the  collisionless 
Rayleigh-Taylor  mechanism  using  techniques  in  two  dimensions  that 
Chaturvedi  and  Kaw  (1975  a,b)  used  in  one  dimension.  Ott  [1977]  has 
described  the  rise  of  equatorial  Spread  F bubbles,  in  both  the  colli- 
sional  and  collisionless  Rayleigh-Taylor  regimes,  by  noting  the  analogy 
with  bubbles  in  fluids.  Szuszczewicz  [1977]  has  discussed  the  chemistry 
and  transport  of  equatorial  Spread  F ionospheric  holes. 

It  has  been  pointed  out  (Scannapieco  and  Ossakow,  1976;  and 
Scannapieco  et  al . , 1976)  that  there  is  a considerable  analogy  between 
equatorial  Spread  F (in  the  collision-dominated  regime)  and  barium  cloud 
studies.  Indeed  barium  clouds  can  be  viewed  as  Pedersen  conductivity  or 
plasma  density  enhancements  while  natural  equatorial  Spread  F bubbles 
(or  artificially  created  ionospheric  holes,  e.g.,  see  Mendi  1 lo  et  al . , 

Note:  Manuscript  submitted  November  7,  1977. 
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1975;  and  Bernhardt  and  da  Rosa,  1977)  can  be  viewed  as  Pedersen  con- 


ductivity or  plasma  density  depletions.  In  the  present  paper,  we  wish 
to  discuss  equatorial  Spread  F bubble  shapes  and  vertical  rise  rates 

within  the  context  of  the  collisional  Rayleigh-Taylor  regime  (as  invoked 
by  Scannapieco  and  Ossakow,  1976),  making  use  of  the  electrodynamic 
analogy  (i.e.,  the  electrical  properties  of  the  bubbles)  with  barium 
clouds.  In  section  II  we  present  the  general  theory,  with  subsections 
dwelling  on  linear  and  nonlinear  bubble  rise  velocity  models.  Section 
III  deals  with  the  discussion  and  summary. 

II.  THEORY 


The  equations  governing  equatorial  Spread  F (within  the  context  of 
the  collisional  Rayleigh-Taylor  instability)  are  (see  Scannapieco  and 
Ossakow,  1976) 

3n 


3t  * 


V t|;  X b) 
1 


= -VR(n-n^) 


(1) 


V^(VinnV^4))  = | a.  ' + |(a.  x b)  • V^n  (2) 

where  c is  the  speed  of  light,  B the  magnetic  field  (taken  to  be  con- 
stant), e the  electronic  charge,  M the  mass  of  the  ions,  £ acceleration 

of  gravity,  b = B/|B|  (pointing  North),  the  subscript  i means  perpendi- 
^ / 
cular  to  b (Note;  the  geometry  we  will  be  employing  for  equatorial 

Spread  F is  such  that  x is  upward,  y is  horizontal  and  eastward  and  z 
is  northward  and  along  the  earth's  magnetic  field),  Vj^  is  the  recombina- 
tion rate,  n the  ion  number  density  (or  equivalently  a field  line  inte- 
grated number  density  in  a slab  geometry),  E = -ViJj,  where  4)  is  the 
electrostatic  potential  (in  this  model  there  is  no  ambient,  i.e.,  zero 
order,  horizontal  electric  field),  the  ion-neutral  collision 
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frequency  (a  function  of  altitude)  and  n^  the  equilibrium  ionospheric 
electron  (ion)  number  density. 

If  we  neglect  recombination  effects,  neglect  the  variation  of 
with  altitude,  and  set  = 4)^  + cf)  (where  from  (1)  and  (2)  Vc})^  = ^) , 
then  we  have 


|i-  X b)  • V.  n = 0 (3) 

V • (nV  4))  = ^ g X b • V n (4) 

1 1 ' e -2-  1 ' ' 

(Note:  The  same  potential  equation  (4)  can  be  obtained  by  keeping  terms 
of  o(l)  compared  with  where  = eB/Mc,  in  eq.  (2)).  The  simplest 

set  of  equations  describing  F region  barium  cloud  morphology  and  stria- 
tions  (for  a one  level  model  and  in  the  x B drift  frame)  is  (see 
Perkins  et  al . , 1973  and  Zabusky  et  al . , 1973) 


xb)  ■ V_N  = 0 


(5) 


V (N  V 4))  = E • V N 
1 ' 1^'  —0  1 


(6) 


where  is  taken  to  be  the  ambient  electric  field  perpendicular  to  B 
(£  E - v^4))  and  N is  the  magnetic  field  line  integrated  Pedersen 
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conductivity  or  equivalently  the  magnetic  field  line  integrated  plasma 
density.  (Note;  In  these  equations  is  taken  to  be  a constant  and 
so  represents  a rotated  geometry,  i.e.,  the  barium  cloud  equations  are 
more  applicable  to  high  latitudes  where  B is  vertical,  V N and  E are 

— X —0 

horizontal.)  If  gravitational  effects  were  also  included  in  (6)  then 
E;*  = (M/e)  £ x b would  replace  E^  on  the  right  hand  side  of 

(6)  (see  Perkins  et  al . , 1973).  The  two  sets  of  equations  (3)  and  (4), 
and  (5)  and  (6)  are  essentially  the  same  with  (fl/v^^)  (M/e)  g^  x b playing 
the  role  of  an  electric  field  in  the  potential  eg.  (4)  (recall  also  the 
definition  of  E^*) . 

Let  us  now  discuss  some  depletion  and  enhancement  cloud  morphology. 
The  separate  coupled  sets  of  equations  (3)  and  (4)  and  (5)  and  (6)  have 
no  steady  state  solutions  if  the  cloud,  as  described  by  n or  N,  has  a 
finite  size  (with  the  gradient  of  n or  N finite)  in  two  dimensions  per- 
pendicular (x,y)  to  the  magnetic  field  (see  Perkins  et  al . , 1973  and 
Dungey,  1958).  Another  result  concerns  the  initial  deformation  of  a 
cylindrical  plasma  density  or  Pedersen  conductivity  enhancement  (barium 
cloud)  and  a cylindrical  plasma  density  or  Pedersen  conductivity  deple- 
tion (equatorial  Spread  F bubbles  or  artificial  ly  created  ionospheric 
holes).  The  situation  is  illustrated  in  Figure  1.  Here  in  this  context 
eq.  (4)  or  (6)  is  simply  the  problem  of  cylindrical  dielectric  immersed 

in  a uniform  electric  field  (where  in  the  equatorial  Spread  F situation 
/\ 

(J2/v^.^) (M/e)  £ X b plays  the  role  of  an  equivalent  electric  field  . In 
the  case  (A)  of  the  density  enhancement  (e.g.,  barium  cloud)  the  induced 
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electric  field  is  opposed  to  the  uniform  electric  field,  E^,  so  that  , 

the  total  electric  field  inside  the  enhancement  is  less  than  E . In  the 

-0 

case  (B)  of  the  density  depletion  (e.g.,  equatorial  Spread  F bubble  or 
artificially  created  "hole")  the  induced  electric  field  is  in  the  direc- 
tion of  the  uniform  electric  field,  E^,  so  that  £' , the  total  electric 
field  inside  the  depletion  is  greater  than  (Note:  the  formulas  for 

are  approximate  and  specific  models  will  be  discussed  later  on  in  the 
text.)  It  should  be  pointed  out  that  in  Fig.  la  the  circles  representing 
the  isodensity  contours  are  such  that  the  enhancement  and  depletion  are 
m.aximum  in  the  center  and  fall  off  as  one  moves  away  from  the  center. 

In  this  way  the  induced  electric  field  goes  to  zero  outside  of  the 
enhancement  or  depletion  (Note:  for  the  enhancement  the  induced  electric 
field  is  — 6N^/(N^  + 6N),  where  N^  is  the  ambient  ionospheric  density 
and  6N  is  the  enhanced  density,  while  for  the  depletion  the  induced 
electric  field  is  6nE^/(n^  - 6n)  where  6n  is  the  magnitude  of  the  deple- 
tion and  n^  is  the  ambient  density).  In  case  (A)  the  x drift  of 
the  center  of  the  enhancement  counteracts  the  general  x ^ drift  and 
the  center  of  the  enhancement  (cloud)  has  the  slowest  drift.  The  back- 
side (side  opposite  to  E^  x ^motion)  of  the  enhancement  catches  up  with 
the  center  and  a steep  density  gradient  forms  there  (depicted  in  Fig.  lb). 
Similarly,  the  frontside  elongates.  The  area  enclosed  by  any  isodensity 
contour  remains  constant  during  the  motion  and  the  maximum  value  of  |N| 
is  a constant  of  the  motion  (incompressibility). 

In  case  (B),  the  x ^ drift  of  the  center  of  the  depletion  (bubble) 
acts  in  concert  with  the  general  x ^ drift  and  the  center  of  the 
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depletion  has  the  fastest  drift.  The  center  now  catches  up  to  the  front- 
side and  a steep  density  gradient  forms  there  (depicted  in  Fig.  lb).  In 
this  case  the  so-called  backside  (using  barium  cloud  terminology)  elongates. 
In  the  context  of  eq.  (3)  and  (4)  in  Fig.  1,  the  area  enclosed  by  any 
isodensity  contour  remains  constant  during  the  motion  and  the  maximum 
value  of  |n|  is  a constant  (incompressibility).  In  both  the  enhancement 
and  depletion  cases  the  density  piles  up  in  the  direction  formed  by  the 
induced  electric  field  crossed  with  the  magnetic  field.  In  both  cases 
where  these  steep  gradients  form  instabilities  should  be  generated.  In 
some  sense  the  density  depletion  can  be  likened  to  an  "anti barium  cloud." 
Within  the  context  of  Fig.  1 we  have  spoken  about  depletions  in  general 
and  the  general  results  apply  to  natural  Spread  F bubbles  or  artificially 
created  ionospheric  holes  (using  H2O,  H2.  or  high  explosives).  Let  us 
now  specialize  these  results  to  equatorial  Spread  F bubbles.  In  this 
case,  neglecting  for  the  moment  any  horizontal  ambient  electric  field, 
in  Fig.  1 is  really  (f^M/ev^^)  £x  b and  it  is  the  polarization 
(induced)  electric  field  (in  the  direction  of  £'  and  equal  to  (fiM/ev^^) 
(<Sn/(nQ  - 6n))  £x  b which  causes  the  bubble  to  rise  with  respect  to 
the  ionosphere  via  the  local  convective  drift  motion.  The  picture 

exhibited  in  Fig.  1 for  Spread  F bubbles  (using  essentially  eq.  (4) 
which  contains  collisional  Rayleigh-Taylor  effects)  shows  that  the 
bubbles  spread  more  vertically  than  horizontally.  This  agrees  with  the 
radar  observations  [Woodman  and  La  Hoz , 1976]  and  of  course  the  numerical 
simulations  [Scannapieco  and  Ossakow,  1976]  which  utilized  eq.  (1)  and  (2). 
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Linear  Spread  F Bubble  Rise  Velocity  Model 
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Let  us  now  look  at  Spread  F bubble  rise  velocities  for  various 
bubble  models.  We  will  utilize  eq.  (4)  and  in  our  first  case  let  us 
include  a horizontal  ambient  electric  field,  E^,  to  show  its  effects. 
Then  the  equation  becomes 


V (nV  (j))  = E*  • V n 

i'  1^'  1 


(7) 


fi  M 

- £ -|>  LL  q ^ 

- -0  e 


Our  first  model  will  be  a linearized  model  such  that 


n = n^(x)  + nie^ky-wt) 


(8) 


(j)  - (j) 


1® 


i (ky-wt) 


where  n^(x)  is  the  equilibrium  altitude  dependent  density  profile  and 
subscript  one  refers  to  perturbed  quantities.  Substituting  (8)  into  (7) 
and  taking  E^  to  be  in  the  y direction  yields  for  the  polarization 
potential 


(j),  = -i 


i + E 

l^in  " 


\ 

o)  k n^ 


The  bubble  velocity  is 


(9) 
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(10) 


In  eq.  (10)  we  use  £=  -Vcjjj  = - so  that  for  the  linear  case  |Vg|  = 
|cE/B|  = |ck(t)j/B|  and  we  have  for  the  vertical  velocity  (6n  e n^) 


(11) 


(Note:  Here  6n  is  positive  and  taken  to  mean  the  magnitude  of  the 
depletion.)  It  is  interesting  to  note  that  adding  in  the  direction 


parallel  (antiparallel)  to  £ x b increases  (decreases)  the  polarization 
field  (as  given  by  -Vcfjj)  and  makes  the  bubbles  rise  faster  (slower).  The 
numerical  simulations  of  Scannapieco  and  Ossakow  [1976]  did  not  have  a 
horizontal  electric  field  and  adding  one  in  the  proper  direction  would 
make  the  bubble  rise  faster.  However  let  us  look  at  the  effect  of  E^. 

We  wish  to  compare  the  two  terms  inside  the  parenthesis  of  eq.  (11).  If 
we  assume  E^  'v-  ImV/m,  B = 0.3  gauss  then  cE^/B  'v-  33  m/sec.  At  350  km 
altitude  g/v^^  40  m/sec,  whereas  at  450  km  g/v^^  'v-  220  m/sec.  Thus 

we  see  that  as  one  goes  to  higher  altitudes  E^  makes  a more  negligible 
contribution  to  the  bubble  rise  velocity.  Indeed,  even  if  the  electric 
field  reverses,  once  the  bubble  gets  to  higher  altitudes  (decreasing 
Vin)  the  bubble  would  keep  rising.  This  bubble  rising  even  after  the 
electric  field  reverses  has  been  seen  experimentally  [Kelley,  private 
communication,  1976]. 
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(12) 


Neglecting  in  eq.  (11)  we  have 


V = _5_  ^ 
B V.  n„ 
in  0 


Using  Fig.  2 for  as  a function  of  altitude,  Fig.  3 depicts  Vg  as  a 
function  of  altitude  for  various  values  of  percentage  depletions,  Sn/n^ 
(large  <Sn/n^  represents  usage  of  eq.  (12)  beyond  the  range  of  validity 
and  perhaps  represents  a lower  bound).  Inmediately  it  can  be  seen  from 
this  figure  that  large  bubble  velocities  ('^^lOO  m/sec)  requires  large 
depletions  and/or  high  altitudes. 

Nonlinear  Spread  F Bubble  Rise  Velocity  Model 

In  the  case  of  a so-called  sheet  geometry  in  which  we  allow  n and 
(p  to  be  functions  of  y only,  eq.  (4)  can  be  integrated  once  to  yield 


3())  _ 

9y  " 


n 


n K 


V 


in 


e ^ 


(13) 


where  n is  the  ambient  density  outside  the  bubble.  For  this  case  the 
0 

dimension  of  the  bubble  in  the  vertical  direction  is  extremely  long 
compared  with  its  horizontal  dimension.  Using  eq.  (13),  the  vertical 
bubble  rise  velocity  now  becomes 


6n 

n 


, 6n 
n 


V. 

in 


(14) 
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where  we  have  used  Vg  = y x b and  n = - 6n  for  a depletion.  A 

peculiarity  of  this  nwdel  is  that  Vg  as  fin/n^  1,  i.e.;  for  a 100% 

depletion.  Fig.  4 shows  the  bubble  rise  velocity  as  a function  of  alti- 
tude for  various  vales  of  6n/n^,  using  eq.  (14)  and  values  of  depicted 
in  Fig.  2.  Here,  as  in  Fig.  3,  high  altitudes  and/or  large  depletions 
yield  high  vertical  bubble  velocities. 

As  we  have  stated  previously,  eq.  (3)  and  (4)  in  general  do  not 
admit  of  a two  dimensional  (i  B^)  solution.  However,  there  is  a patho- 
logical case  where  one  can  obtain  exact  solutions  and  we  will  discuss 
this  now.  We  write  for  the  density 


n(x,y)  = n„  - nn(x,y)  = n. 


1 - ^F(x.y) 
0 


(15) 


where  n^  is  the  ambient  ionospheric  density,  n^  represents  the  depletion 
or  bubble  (fin/n^  is  the  magnitude  of  the  ratio  of  maximum  depleted  density 
to  ambient  density),  and  F describes  the  shape  of  depletion.  If  ng(x,y) 

0 as  x,y  “ and  n is  finite  then  there  is  no  equilibrium  solution 
to  eq.  (3)  and  (4).  However,  if  the  depletion  density  is  a constant  and 
the  shape  of  the  density  depletion  contour  is  an  ellipse  described  by 


ng  = 6n  H 


1 - ^ 


[where  H(x)  is  the  Heaviside  step  function  and  a and  b represent  the 
vertical  dimension  (axis  of  the  ellipse)  and  horizontal  (east-west) 
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dimension  respectively],  then  the  electric  field  inside  the  depletion 
is  a constant  and  the  depletion  (bubble)  moves  with  constant  velocity 
without  changing  its  shape.  This  method  was  first  employed  by  Linson 
[1972]  for  barium  clouds  (plasma  density  enhancements)  where  eq.  (5)  and 
(6)  were  utilized.  (Note:  In  this  case  the  minus  signs  in  eq.  (15) 
would  become  plus  signs.)  Indeed,  eq.  (4)  or  (6)  can  be  thought  of  as 
the  case  of  a dielectric  immersed  in  a uniform  electric  field.  However, 
the  plasma  depletion  is  analogous  to  the  case  of  a cavity  immersed  in  a 
dielectric  with  a uniform  field  (since  the  field  is  enhanced  inside  a 
depletion).  The  solution  of  the  potential  eq.  (4)  or  (6)  for  the  general 
elliptical  shape  can  be  found  in  Smythe  [1950].  For  the  case  of  a plasma 


depletion  described  by  eq.  (4),  with  x b in  the  y direction,  we  have 
for  the  polarization  (induced)  electric  field 


(16) 


with  the  attendant  bubble  rise  velocity  given  by 


(17) 


11 


(Note:  In  arriving  at  eq.  (16)  we  have  assumed  to  be  a constant. 
This  is  a good  approximation  so  long  as  the  size  of  the  depletion  is 
small  compared  with  an  atmospheric  scale  height). 

Now  we  will  examine  eq.  (17)  with  respect  to  different  shapes.  For 
the  case  b « a,  i.e.,  the  sheet  geometry,  we  obtain 


Vn  = 


_9_ 

^in 


§n 

n. 


I M 
’ % 


which  is  eq.  (14)  as  it  should  be.  For  a/b  ->-0,  i.e.,  a slab  case 

/\ 

(uniform  in  the  direction  of  £ x b) , Vg  ->  0 since  the  induced  electric 
field  goes  to  zero.  When  a = b (cylindrical  case),  we  have 


6n 

n„ 


Vo  = 


2 u 1 1 


§r\ 

n. 


_9_ 

^in 


(18) 


For  Sn/n^  « 1 eq.  (18)  yields  Vg  = 1/2  (g/v^^) (6n/n^)  which  is  half 
the  value  given  by  the  linear  perturbation  model  (see  eq.  (12)).  With 
6n/nQ  = 1 (100%  depletion),  eq.  (18)  yields  Vg  = FiO-  5 shows 

the  vertical  bubble  rise  velocity  for  a cylindrical  shape  (as  given  by 
eq.  (18))  and  Fig.  6 exhibits  the  bubble  rise  velocity  for  an  elliptical 
shape  (a/b  = 5)  as  a function  of  altitude  for  various  values  of  6n/n^, 
using  values  of  depicted  in  Fig.  2. 
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III.  DISCUSSION  AND  SUMMARY 


Figures  3-6  depict  vertical  bubble  (plasma  depletion)  rise  veloci- 
ties for  various  bubble  models  (shapes),  for  equatorial  Spread  F geo- 
metry, as  a function  of  altitude  and  percentage  depletions.  Exceedingly 
high  velocities  can  be  attained,  especially  for  the  severely  elongated 
models  (e.g.,  sheet  model  in  Fig.  4 and  the  elliptical  model  with  a 5:1 
axial  ratio  in  Fig.  6).  Table  1 shows  the  bubble  rise  velocity,  in 
units  of  g/v^^,  as  a function  of  fractional  depleted  plasma  density  for 
various  bubble  shapes.  For  all  values  of  Sn/n^  (except  at  1.0)  the  non- 
linear cylindrical  model  gives  a smaller  velocity  (see  also  eqs.  (12) 
and  (18))  than  the  linear  model.  The  table  also  shows  that  the  cylin- 
drical bubble  velocity  for  a given  <Sn/n^  is  attained  by  the  sheet  model 
for  a value  of  1/2  Sn/n^  (see  also  eqs.  (14)  and  (18)).  An  elliptical 
model  with  10:1  axial  ratio  is  also  shown  in  Table  1.  Except  for 
6n/nQ  = 1,  it  has  already  become  fairly  sheet-like.  A common  feature  to 
all  models  is  that  the  higher  the  altitude  and/or  the  higher  the  density 
depletion,  the  higher  the  vertical  rise  velocity.  This  result  is  con- 
sistent with  the  recent  satellite  in  situ  measurements  described  by 
McClure  et  al . [1977]. 

The  models  described  by  Fig.  3-6  essentially  assume  that  the  bubble 
does  not  change  its  shape.  For  real  bubbles,  where  the  depletion  den- 
sity changes  inside  the  bubble,  the  bubble  should  elongate  in  the 
vertical  direction  (steepening  on  the  topside  of  the  bubble)  and  become 
narrower  in  the  east-west  extent.  This  is  in  agreement  with  the  radar 
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observations  (Woodman  and  La  Hoz,  1976)  and  the  numerical  simulation  of 
rising  equatorial  Spread  F bubbles  (Scannapieco  and  Ossakow,  1976). 
Therefore,  in  modeling  a changing  bubble  by  one  of  constant  shape  we  are 
employing  somewhat  of  an  artifact.  However,  using  this  technique  Linson 
[1972]  has  shown  (using  optical  data)  that  evolving  barium  clouds  (plasma 
density  enhancements)  tend  to  be  more  cylindrical  than  sheet-like.  Using 
the  constant  shape  approach  and  the  simulation  results  of  Scannapieco 
and  Ossakow  [1976],  where  Vg  10  m/sec  and  Sn/n^  40%  for  an  altitude 
of  350  km,  the  cylindrical  bubble  shape  is  favored  (although  the  simula- 
tion pictures  show  an  elongated  bubble).  More  data  and  simulation 
results  are  needed  to  ascertain  which  model  fits  best. 

All  of  the  models  described  in  this  paper,  based  on  collisional 
Rayleigh-Taylor  type  equations,  show  that  the  vertical  bubble  rise 
velocity  is  of  the  form 


(19) 


where  f(6n/nQ)  is  some  function  of  6n/n^.  The  factor  9/^^^-^  should  not 
be  too  surprising  since  the  growth  rate  (y)  for  the  collisional  Rayleigh- 
Taylor  instability,  neglecting  recombination  effects,  is  g/v^-^  L (Hudson 
and  Kennel , 1975;  Scannapieco  and  Ossakow,  1976),  which  is  a velocity 
divided  by  the  electron  density  gradient  scale  length.  This  of  course 
is  similar  to  the  growth  rate  (y  = cE^/BL)  for  barium  cloud  striations 


14 


[ 


due  to  X ^gradient  drift  instability  (Simon , 1963;  Linson  and  Workman, 
1970).  This  similarity  between  the  gravitational  and  the  electric  field 
terms  has  been  shown  in  eq.  (11).  In  eq.  (19),  the  factor  g/^^p  favors 
high  altitude  because  decreases  as  a function  of  altitude  (see 
Fig.  2).  The  factor  f increases  with  increasing  6n/n^  (see  Fig.  3-6). 

At  least  according  to  linear  theory,  <5n/n^  increases  with  altitude  and 
decreasing  bottomside  electron  density  gradient  scale  length  (y  = g/v^^  L). 
These  ideas  indicate  that  there  is  a range  of  bubble  vertical  rise 
velocities  (in  accord  with  HcCl ure  et  al . , [1977]),  but  these  rise 
velocities  depend  on  ambient  equatorial  F region  ionospheric  conditions 
(e.g.,  height  of  the  F peak  and  bottomside  gradient  scale  length). 

Recent  numerical  simulations  (S.  Ossakow,  S.  Zalesak  and  B.  E.  McDonald, 
private  communication,  1977)  support  these  ideas.  Also,  it  should  be 
pointed  out  that  for  <5n/n^  = 0,  i.e.,  for  bubbles  which  have  the  same 
density  as  the  background  ionosphere,  the  induced  (polarization)  electric 
field  (which  causes  the  bubble  to  rise  through  £ x ^)  becomes  zero  (with 
all  models,  as  it  should),  as  does  the  attendant  vertical  bubble  rise 
velocity.  This  simply  states  that  a bubble  created  on  the  bottomside  of 
the  F region  will  rise  until  its  density  equals  the  topside  background 
ionospheric  density. 

It  should  be  pointed  out  that  our  vertical  velocity  result  for 
cylindrical  bubbles  does  not  agree  with  Qtt's  [1977]  result.  In  his 
analogy  with  bubbles  in  a fluid,  Ott  found  only  a cylindrical  solution. 

For  the  collision  dominated  regime  he  found  Vg  = l/2(g/v^.^) (6n/n^) , in 
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our  terminology;  whereas,  our  result  shows  Vg  = l/2(g/v^.^) (dn/n^)  • 

(1  - l/2(6n/n^))  For  the  case  where  the  bubble  density  is  negligibly 
small  (6n/n^=:  1)  compared  with  the  density  in  the  ionosphere  outside, 
Ott's  result  yields  Vg  = l/2(g/v^^);  whereas,  our  result  yields 

''b  = 

Finally  some  comments  about  the  various  upward  and  westward  bubble 
velocities  presented  by  McCl ure  et  al . [1977]  and  the  tilts  of  the 
plumes  seen  by  the  Jicamarca  radar  measurements  (Woodman  and  La  Hoz, 

1976).  Equations  (1)  and  (2),  used  by  Scannapieco  and  Ossakow  [1976], 
or  eqs.  (3)  and  (4)  for  that  matter  produce  only  vertical  bubble  rise. 
These  equations  only  admit  of  a horizontal  polarization  electric  field 
(eastward)  caused  by  the  £ x b term  such  that  depletions  rise  and  enhance- 
ments fall  with  respect  to  the  ionosphere.  Inclusion  of  an  ambient 
background  horizontal  electric  field  (eastward  or  westward)  in  the  equa- 
tions could  only  cause  additional  vertical  (upward  or  downward)  motion 
of  the  bubbles.  Since  the  ionosphere  moves  mainly  downward  at  night 
(v^estward  electric  field;  see  McClure  et  al . , 1977)  this  could  slow  down 
the  bubbles  at  low  altitudes  and  have  a negligible  effect  at  higher 
altitudes  (see  eq.  (11)).  Vertically  the  bubbles  move  more  slowly  at 
lower  altitudes  and/or  small  values  of  6n/n^.  Thus  there  exist  conditions 
whereby  the  bubbles  do  not  move  vertically  with  respect  to  the  iono- 
sphere. McCl ure  et  al . [1977]  and  Woodman  and  La  Hoz  [1976]  have 
stated  that  the  ionosphere  moves  mainly  eastward  at  night.  This  would 
mean  an  ambient  electric  field  pointing  downward.  Inclusion  of  such  an 
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electric  field  in  the  equations  would  cause  bubbles  to  move  eastward. 

Thus  the  westward  motion  of  the  bubbles  must  be  caused  by  introducing 
an  asymmetry  into  the  eq.  (1)  and  (2)  (or  (3)  and  (4))  which  causes  a 
vertical  (upward)  electric  field.  Such  an  effect  can  be  produced  by  an 
eastward  directed  neutral  wind  which  would  essentially  polarize  the 
bubble  to  produce  an  upward  polarization  electric  field  and  account  for 
the  westward  motion  of  the  bubble  (and  tilts  of  the  radar  plumes).  This 
last  idea  has  also  been  set  forth  by  V/oodman  and  La  Hoz  [1976]  and  Ott 
[1977]. 

Recalling  our  discussion  about  the  electric  field,  E*  h E„  + (il/v.  ). 

3 ’-o'  in' 

(M/e)  £ X b,  after  eq.  (6),  all  of  our  vertical  bubble  rise  velocities  can 
take  account  of  horizontal  motions.  The  general  bubble  velocity  is  then 
just  Vg  = (c/B)  X b f(6n/nQ)  = ((c/B)  ^ x b - a/v^^)  f(6n/nQ).  This 
clearly  shows  how  the  directionality  of  Vg  is  formed.  In  this  case  for 
the  elliptical  shape  a is  in  the  direction  of  x b,  etc.  and  non- 
constant density  depletions  will  steepen  and  elongate  in  this  direction. 

Recent  results  from  the  Wideband  associated  equatorial  experiments 
performed  at  Kwajalein  (August  1977)  lend  credence  to  this  notion  (S . 
Ossakow  and  J.  Clark,  private  communication,  1977).  Backscatter  irre- 
gularity ('^'l  meter)  measurements  (performed  by  D.  Towle)  have  shown  that 
these  patches  rise  and  appear  westward  when  neutral  wind  measurements 
(Fabry -Perot  interferometer  measurements  performed  by  M.  Biondi  and  D.  Sipler)  show 
strong  west  to  east  winds  ('^^160  m/sec)  earlier  in  the  evening.  As  the 
evening  progresses  the  neutral  winds  die  down  ('^^10-20  m/sec)  and  the 
radar  shows  the  irregularities  appearing  eastward. 
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Table  1.  Bubble  rise  velocity  (in  units  of  Vg,  as  a function 

of  fractional  depleted  density,  6n/nQ,  for  various  bubble  shapes. 
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Fig.  1 — Deformation  of  (A)  a plasma  density  enhancement  (e.g.,  a barium 
cloud)  and  (B)  a plasma  density  depletion  (e.g.,  equatorial  Spread  F bubble 
or  artificially  created  “hole”)  in  crossed  electric  (^)  and  magnetic  (B)  fields 
(representative  of  equatorisJ  geometry),  (a)  Initial  (t=0)  cylindrical  isoden- 
sity contours  with  largest  enhancement  (+)  and  depletion  (-)  values  in  the 
center  and  falling  off  away  from  these  values  to  an  ambient  level  n^  outside. 

represents  the  total  electric  field  inside  the  enhancement  or  depletion 
and  ^ is  the  ambient  field  outside,  (b)  Qualitative  evolution  of  the  isoden- 
sity contours,  (c)  Approximate  values  for  the  initial  electric  field  (E')  inside 
the  enhancement  or  depletion. 


23 


Fig.  2 — Ion-neutral  collision  frequency,  as  a 
function  of  altitude 
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Fig.  3 — Depletion  (bubble)  vertical  rise  velocity,  Vg, 
as  a function  of  altitude  for  various  values  of  the  per- 
centage depletion,  Sn/n^,  for  the  linear  bubble  model, 
i.e.,  Vg  = (glUiJ(8nlnJ 
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ALTITUDE  (km) 


Fig.  4 Same  as  fig.  3,  but  for  the  nonlinear  sheet  bubble 
model,  i.e.,  Vg  = (g/yin)(6n/no)(l-6n/no)-l 
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ELLIPTICAL  M0DEL(5:i)BUBBLE  VEOXITY.  (m/sec) 


Fig.  6 — Same  as  fig.  3,  but  for  the  nonlinear  elliptical 
bubble  model  (with  the  vertical  dimension  five  times  the 
horizontal  dimension),  i.e.,  Vg  = (g/u^)(5  5n/nQ)(6-5 
5n/no)-i 
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OICY  ATTN  DRSEL-BL-SY-S  F.  E.  NILES 


DIRECTOR 

BMD  ADVANCED  TECH  CTR 
HUNTSVILLE  OFFICE 
P.  0.  BOX  1500 
HUNTSVILLE,  AL  35807 
OICY  ATTN  ATC-T 
OICY  ATTN  ATC-0 
OICY  ATTN  ATC-R 


MELVIN  T.  CAPPS 
W.  DAVIES 
DON  RUSS 


PROGRAM  MANAGER 
BMD  PROGRAM  OFFICE 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

OICY  ATTN  DACS-BMT  JOHN  SHEA 


CHIEF  C-E  SERVICES  DIVISION 
U.S.  ARMY  COMMUNICATIONS  CMD 
PENTAGON  RM  1B269 
WASHINGTON,  D.C.  20310 

OICY  ATTN  CC-OPS-CE 
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COMMANDER 

HARRY  DIAMOND  LABORATORIES 
2800  POWDER  MILL  ROAD 
ADELPHI,  MD  20783 

(CNWDI-INNER  ENVELOPE:  ATTN:  DRXDO-RBH) 
OICY  ATTN  MILDRED  H.  WEINER  DRXDO-TI 
OICY  ATTN  DRXDO-RB  ROBERT  WILLIAMS 
OICY  ATTN  DRXDO-NP  FRANCIS  N.  WIMENITZ 
OICY  ATTN  DRXDO-NP  CYRUS  MOAZED 


DIRECTOR 

TRASANA 

WHITE  SANDS  MISSILE  RANGE,  NM  88002 
OICY  ATTN  ATAA-SA 
OICY  ATTN  TCC/F,  PAYAN  JR 
OICY  ATTN  ATAA-TAC  LTC  JOHN  HESSE 


COMMANDER 

U.S,  ARMY  COMM-ELEC  ENGRG  INSTAL  AGY 
FT.  HUACHUCA,  AZ  85613 

OICY  ATTN  EED-PED  GEORGE  LANE 


COMMANDER 

U.S.  ARMY  ELECTRONICS  COMMAND 
FORT  MONMOUTH,  NJ  07703 

OICY  ATTN  DRSEL-NL-RD  H.  S.  BENNET 
OICY  ATTN  DRSEL-PL-ENV  HANS  A.  BOMKE 


COMMANDER 

U.S.  ARMY  FOREIGN  SCIENCE  & TECH  CTR 
220  7TH  STREET,  NE 
CHARLOTTESVILLE,  VA  22901 

OICY  ATTN  P.  A.  CROWLEY 
OICY  ATTN  R.  JONES 


COMMANDER 

U.S.  ARMY  MATERIEL  DEV 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 
OICY  ATTN  DRCLDC 


& READINESS  CMD 


J.  A.  BENDER 
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COMMANDER 

U.S.  ARMY  NUCLEAR  AGENCY 
FORT  BLISS,  TX  79916 

OICY  ATTN  MONA-WE  J.  BERBERET 


DIRECTOR 

U.S.  ARMY  BALLISTIC  RESEARCH  LABS 
ABERDEEN  PROVING  GROUND,  MD  21005 
OICY  ATTN  LAWRENCE  J.  PUCKETT 


COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONMOUTH,  NJ  07703 

OICY  ATTN  DOCUMENT  CONTROL 


COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35809 
OICY  ATTN  JIM  GAMBLE 


CHIEF  OF  NAVAL  OPERATIONS 
NAVY  DEPARTMENT 
WASHINGTON,  D.C.  20350 

OICY  ATTN  OP  943  LCDR  HUFF 
OICY  ATTN  ALEXANDER  BRANDT 
OICY  ATTN  RONALD  E.  PITKIN 


CHIEF  OF  NAVAL  RESEARCH 
NAVY  DEPARTMENT 
ARLINGTON,  VA  22217 

OICY  ATTN  CODE  418 
OICY  ATTN  CODE  461 


COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
NAVAL  ELECTRONIC  SYSTEMS  CMD  HQS 
WASHINGTON,  D.C.  20360 

OICY  ATTN  NAVALEX  034  T BARRY  HUGHES 

OICY  ATTN  PME  117-T  SATELLITE  COMM  PROJECT  OFF 

OICY  ATTN  PME  117 


COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CTR 
4301  SUITLAND  ROAD,  BLDG.  5 
WASHINGTON,  D.C.  20390 

OICY  ATTN  MR.  DUBBIN  STIC  12 
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COMMANDER 


NAVAL  ocean  systems  CENTER 

SAN  DIEGO, 

CA  92152 

0 ICY 

ATTN  WILLIAM  F . 

MOLER 

0 ICY 

ATTN  CODE  0230 

C.  BAGGETT 

03CY 

ATTN  CODE  2200 

0 ICY 

ATTN  R.  EASTMAN 

DIRECTOR 

NAVAL  RESE 

ARCH  LABORATORY 

WASHINGTON 

, D.C.  20375 

OICY 

ATTN  CODE  7700 

TIMOTHY  P.  COFFEY 

03CY 

ATTN  CODE  7701 

JACK  D.  BROWN 

0 ICY 

ATTN  CODE  7750 

OICY 

ATTN  CODE  5400 

BRUCE  WALD 

03CY 

ATTN  CODE  5450 

TELECOMMUNICATIONS  SYSTEMS 

COMMANDER 

NAVAL  SPAC 

E SURVEILLANCE  S 

YSTEM 

DAHLGREN, 

VA  22448 

OICY 

ATTN  CAPT  J.  H. 

BURTON 

COMMANDER 

NAVAL  SURF 

ACE  WEAPONS  CENTER 

WHITE  OAK, 

SILVER  SPRING, 

MD  20910 

0 ICY 

ATTN  CODE  WA501 

NAVY  NUC  PRGMS,  OFF 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
NAVY  DEPARTMENT 
WASHINGTON,  D.C.  20376 
OICY  ATTN  NSP-2141 

OICY  ATTN  NSSP-2722  FRED  WIMBERLY 


NAVAL  SPACE  SYSTEM  ACTIVITY 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CALIF.  90009 

OICY  ATTN  A.  B.  HAZZARD 


COMMANDER 

ADC/DC 

ENT  AFB,  CO  80912 

OICY  ATTN  DC  MR.  LONG 


-c- 


COMMANDER 

ADCOM/XPD 

ENT  AF3,  CO  80912 

OICY  ATTN  XPQDQ 


AF  GEOPHYSICS  LABORATORY.  AFSC 


MANSCOM  AFB,  MA  017 
OICY  ATTN  OPR 
OICY  ATTN  OPR 
OICY  ATTN  LKB 
OICY  ATTN  OPR 
OICY  ATTN  SUOL 
OICY  ATTN  PHP 
OICY  ATTN  PHD 
OICY  ATTN  PHD 


1 

HAROLD  GARDNER 
JAMES  C.  ULWICK 
KENNETH  S.  W.  CHAMPION 
ALVA  T.  STAIR 
RSCH  LIB 
JULES  AARONS 
JURGEN  BUCHAU 
JOHN  P.  MULLEN 


AF  WEAPONS  LABORATORY, 
KIRTLAND  AFB,  NM  87117 


AFSC 


AFTAC 


OICY 

ATTN 

SUL 

0 ICY 

ATTN 

CA 

ARTHUR  H.  GUENTHER 

0 ICY 

ATTN 

DYC 

CAPT  L.  WITTWER 

0 ICY 

ATTN 

SAS 

JOHN  M.  KAMM 

0 ICY 

ATTN 

DYC 

CAPT  MARK  A.  FRY 

CK  AFE 

1,  FL 

3 29  25 

0 ICY 

ATTN 

TF/MAJ  WILEY 

0 ICY 

ATTN 

TN 

AIR  FORCE  AVIONICS  LABORATORY,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AFAL  AAB  H.  M.  HARTMAN 
OICY  ATTN  AAD  ALLEN  JOHNSON 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/XR 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  XRC  LTC  J.  MORIN 
OICY  ATTN  XRE  LT  MICHAELS 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/YS 
HANSCOM  AFB,  MA  01731 
OICY  ATTN  YSEV 
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COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  NICD  LIBRARY 
OICY  ATTN  ETD  B.  L.  BALLARD 


HQ  USAF/RD 

WASHINGTON,  D.C.  20330 
OICY  ATTN  RDQ 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13440 

OICY  ATTN  EMTLD  DOC  LIBRARY 
OICY  ATTN  OCSE  V.  COYNE 


SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

OICY  ATTN  SZJ  MAJOR  LAWRENCE  DOAN 


COMMANDER  IN  CHIEF 
STRATEGIC  AIR  COMMAND 
OFFUTT  AFB,  NB  68113 

OICY  ATTN  XPFS  MAJ  BRIAN  G.  STEPHAN 
OICY  ATTN  ADWATE  CAPT  BRUCE  BAUER 
OICY  ATTN  NRT 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION  (AFSC) 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  JIM  DEAS 


SAMSO/ YA 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  YAT  CAPT  L. 


BLACKWELDER 


SAMSO/ 5K 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  SKA  LT  MARIA  A.  CLAVIN 
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SAMSO/MN 

NORTON  AFB,  CA  92419 
CMINUTEMAN) 

OICY  ATTN  MNNL  LTC  KENNEDY 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01733 

OICY  ATTN  ETEI  A.  LORENTZEN 
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U.S.  ENERGY  RSCH  AND  DEV  ADMIN 


EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.  0.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  JAMES  R.  BREEDLOVE 


UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  TECH  INFO  DEPT  L-3 
OICY  ATTN  RONALD  L.  OTT  L-531 
OICY  ATTN  DONALD  R.  DUNN  L-156 
OICY  ATTN  RASPH  S.  HAGER  L-31 


LOS  ALAMOS  SCIENTIFIC  LABORATORY 
P.  0.  BOX  1663 


ALAMOS 

, NM 

87545 

OICY 

ATTN 

DOC 

CON 

FOR 

ERIC 

L INDMAN 

OICY 

ATTN 

DOC 

CON 

FOR 

R.  F 

. TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

ERIC 

JONES 

0 ICY 

ATTN 

DOC 

CON 

FOR 

JOHN 

S.  MALIK 

0 ICY 

ATTN 

DOC 

CON 

FOR 

MARTIN  TIERNEY 

OICY 

ATTN 

DOC 

CON 

FOR 

JOHN 

Z INN 

SANDIA  LABORATORIES 
P.  0.  BOX  5800 


ALBUQUERQUE 

, NM 

87 

115 

OICY 

ATTN 

DOC 

CON 

FOR 

J . 

P . 

MARTIN  ORG  1732 

0 ICY 

ATTN 

DOC 

CON 

FOR 

W. 

D. 

BROWN  ORG  1353 

0 ICY 

ATTN 

DOC 

CON 

■ FOR 

A. 

DEAN  THORNBROUGH  ORG 

OICY 

ATTN 

DOC 

CON 

FOR 

T . 

WRIGHT 

0 ICY 

ATTN 

DOC 

CON 

FOR 

D. 

A. 

DAHLGREN  ORG  1722 

-10- 


OTHER  GOVERNMENT 


DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TELECOMMUNICATIONS 
INSTITUTE  FOR  TELECOM  SCIENCE 
BOULDER,  CO  80302 

OICY  ATTN  WILLIAM  F.  UTLAUT 
OICY  ATTN  G.  REED 


NATIONAL  OCEANIC  6 ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICY  ATTN  JOSEPH  H.  POPE 
OICY  ATTN  RICHARD  GRUBB 
OICY  ATTN  C.  L.  RUFENACH 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 

OICY  ATTN  ATS-6  OFC  P.  CORRIGAN 


-11- 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 


P.  0.  BOX  92957 
LOS  ANGELES,  CA 
OICY  ATTN 

90009 

IRVING 

M.  GARFUNKEL 

OICY 

ATTN 

T . 

M. 

SALMI 

0 ICY 

ATTN 

V . 

JOSEPHSON 

0 ICY 

ATTN 

S . 

P . 

BOWER 

0 ICY 

ATTN 

N . 

D. 

STOCKWELL 

0 ICY 

ATTN 

P . 

P . 

OLSEN  120  RM 

2224E 

0 ICY 

ATTN 

F . 

E . 

BOND  RM  5005 

0 ICY 

ATTN 

J . 

E . 

CARTER  120  RM  2209 

OICY 

ATTN 

F . 

A. 

MORSE  A6  RM 

2407 

ANALYTICAL 

SYSTEMS 

ENGINEERING  CORP 

5 OLD  CONCORD  ROAD 
BURLINGTON,  MA  01805 

OICY  ATTN  RADIO  SCIENCES 


BOEING  COMPANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  WA  98124 

OICY  ATTN  GLEN  KEISTER 
OICY  ATTN  D.  MURRAY 


CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
5175  MIRAMAR  ROAD 
LA  JOLLA,  CA  92037 

OICY  ATTN  HENRY  G.  BOOKER 


BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  55807 

OICY  ATTN  ROMEO  A.  DELIBERIS 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 


555  TECHNOLOGY  SQUARE 

CAMBRIDGE, 

MA  02139 

0 ICY 

ATTN  D.  B. 

COX 

OICY 

ATTN  J.  P. 

GILMORE 

-12- 


COMPUTER  SCIENCES  CORPORATION 
P.  0.  BOX  530 
6565  ARLINGTON  BLVD 
FALLS  CHURCH,  VA  22046 
OICY  ATTN  H.  BLANK 
OICY  ATTN  JOHN  SPOOR 


COMSAT  LABORATORIES 
LI  NTH  I CUM  ROAD 
CLARKSBURG,  MD  20734 

OICY  ATTN  R.  R.  TAUR 


CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.  T.  FARLEY  JR 


ESL  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 

OICY  ATTN  R.  K.  STEVENS 

OICY  ATTN  J.  ROBERTS 

OICY  ATTN  JAMES  MARSHALL 

OICY  ATTN  V.  L.  MOWER 

OICY  ATTN  C.  W.  PRETTIE 


FORD  AEROSPACE  & COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  ALTO,  CA  94303 

OICY  ATTN  J.  T.  MATTINGLEY  MS  X22 


GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 
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GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 


OICY 

ATTN 

DASIAC 

OICY 

ATTN 

DON  CHANDLER 

OICY 

ATTN 

TOM  BARRETT 

0 ICY 

ATTN 

TIM  STEPHENS 

0 ICY 

ATTN 

WARREN  S.  KNAPP 

0 ICY 

ATTN 

WILLIAM  MCNAMERA 

OICY 

ATTN 

B.  GAMBILL 

OICY 

ATTN 

MACK  STANTON 

GENERAL  RESEARCH  CORPORATION 

P.  0.  BOX  3587 

SANTA  BARBARA,  CA  93105 

OICY  ATTN  JOHN  ISE  JR 
OICY  ATTN  JOEL  GARBARINO 


GEOPHYSICAL 

UNIVERSITY 

FAIRBANKS, 

OICY 
OICY 
0 ICY 


INSTITUTE 
OF  ALASKA 
AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFF 
ATTN  T.  N.  DAVIS  (UNCL  ONLY) 
ATTN  NEAL  BROWN  (UNCL  ONLY) 
ATTN  TECHNICAL  LIBRARY 


CER) 


GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHAL  CROSS 


INSTITUTE  FOR 
400  ARMY-NAVY 


DEFENSE 

DRIVE 


ANALYSES 


ARLINGTON, 

VA 

2220  2 

OICY 

ATTN 

J . M . 

AEIN 

OICY 

ATTN 

ERNES 

T BAUER 

OICY 

ATTN 

HANS 

WOLFHARD 

OICY 

ATTN 

JOEL 

BENGSTEN 

HARRIS  CORP,  E. 

S .D. 

P.  0.  BOX 

37 

MELBOURNE, 

FL 

329  0 1 

OICY 

ATTN 

ADV  PROG  DEPT 

DR 


CARL  DAVIS 
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HSS,  INC. 

2 ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 


INTL  TEL  S TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 


JAYCOR 

1401  CAMINO  DEL  MAR 
DEL  MAR,  CA  92014 

OICY  ATTN  S.  R.  GOLDMAN 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 


LOCKHEED  MISSILES  S SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.  R.  CHURCHILL 


LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 

OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 

OICY  ATTN  BILLY  M.  MCCORMAC  DEPT  52-54 


KAMAN  SCIENCES  CORP 

P.  0.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 

OICY  ATTN  B.  J.  BITTNER 


LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 
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M.I.T.  LINCOLN  LABORATORY 
P.  0.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  LIB  A-082  FOR  DAVID  M.  TOWLE 

OICY  ATTN  MR,  WALDEN  X113 

OICY  ATTN  JAMES  H.  PANNELL  L-246 

OICY  ATTN  D.  CLARK 


MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  BILL  OLSON 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


OICY 

ATTN 

P.  FISCHER 

OICY 

ATTN 

W,  F,  CREVIER 

OICY 

ATTN 

STEVEN  L,  GUTSCHE 

OICY 

ATTN 

D,  SAPPENFIELD 

OICY 

ATTN 

R,  BOGUSCH 

OICY 

ATTN 

R,  HENDRICK 

OICY 

ATTN 

RALPH  KILB 

OICY 

ATTN 

DAVE  SOWLE 

OICY 

ATTN 

F.  FAJEN 

OICY 

ATTN 

M,  SCHEIBE 

OICY 

ATTN 

CONRAD  L-  LONGMIRE 

OICY 

ATTN 

WARREN  A,  SCHLUETER 

MITRE  CORPORATION,  THE 
P,  0,  BOX  208 
BEDFORD,  MA  01730 


OICY 

ATTN 

J • 

C.  KEENAN 

OICY 

ATTN 

G. 

HARDI NG 

OICY 

ATTN 

CHIEF  SCIENTIST  W 

OICY 

ATTN 

S , 

A MORIN 

OICY 

ATTN 

C . 

E-  CALLAHAN 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD, 

SANTA  MONICA,  CA  90404 

OICY  ATTN  E,  C.  FIELD  JR 
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PHOTOMETRI CS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

OICY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  21589 
SEATTLE,  VIA  98111 

OICY  ATTN  E.  J.  FREMOUW 


PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

OICY  ATTN  A.  THOMPSON 
OICY  ATTN  JOSEPH  B.  WORKMAN 


R S D ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

OICY  ATTN  RICHARD  LATTER 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT  JR 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 


SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  D.  SACHS 
OICY  ATTN  E.  A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGAL 


RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 
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SCIENCE  APPL /CAT  IONS,  INC. 
HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

OICY  ATTN  DALE  H.  DIVIS 


SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 
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